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A new class of random quantum-dynamical systems in 
continuous space is introduced and studied in some detail. 
Each member of the class is characterized by a Hamiltonian 
which is the sum of two parts. While one part is determinis- 
tic, time-independent, and quadratic, the Weyl-Wigner sym- 
bol of the other part is a homogeneous Gaussian random field 
which is delta correlated in time and arbitrary, but smooth in 
position and momentum. Exact expressions for the time evo- 
lution of both (mixed) states and observables averaged over 
randomness are obtained. The differences between the quan- 
tum and the classical behavior are clearly exhibited. As a 
special case it is shown that, if the deterministic part corre- 
sponds to a particle subjected to a constant magnetic field, 
the spatial variance of the averaged state grows diffusively for 
long times independent of the initial state. 
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The spatial spreading of a state under the free time 
evolution is a well-known and fundamental phenomenon 
in non-relativistic quantum and classical mechanics [Q. 
In order to make a quantitative statement, let of denote 
the variance of the position at time t of a spinless point 
particle moving in continuous space, ft is assumed that 
the particle was prepared initially in some state, which 
is normalized but not necessarily pure. For a free par- 
ticle a simple calculation then shows that of increases 
asymptotically for large t as of ~ t 2 gl/m 2 . Here m > 
is the mass of the particle and Qq is the variance of its 
momentum in the initial state. This relation holds both 
in the quantum and in the classical case. However, clas- 
sical states, e.g. pure ones, may have a sharp momentum, 
that is, Qq = 0, whereas p§ > for all quantum states, 
including the pure ones ( "wave-packet spreading" ) . ft is 
clear that of can also be calculated exactly for a time 
evolution governed by a more general Hamiltonian being 
at most quadratic in momentum and position ||]. 

A variety of physical systems, whose properties are un- 
known in detail, are successfully modeled by a Hamilto- 
nian with a random part added to a (simple) determin- 
istic part. Thus it is a challenging problem to study the 
properties of the random variable of and of related key 
quantities in these models. For recent investigations de- 
voted to the question, how the ballistic long-time behav- 
ior of of of the free particle is modified by adding a Gaus- 



sian time-dependent random potential, see |3|-|5| and ref- 
erences therein. A stepping stone in this direction was 
the work of Jayannavar and Kumar || , who — in building 



on treatments of lattice models Q — exploited the simpli- 
fying feature of a vanishing correlation time ( "Gaussian 
white-noise potential"). Their main result concerns the 
quantum case with a particular pure initial state. They 
derived an exact expression for the spatial variance T, 2 of 
the averaged state at time t and found S| ~ t 3 for large t. 
(Note that the averaged spatial variance of never exceeds 
Sf .) Interesting attempts to incorporate a non-zero cor- 
relation time by the use of perturbative methods — with 
partially conflicting results — can be found in M . 

However, several problems of considerable interest 
have not been tackled so far. Firstly, in order to de- 
scribe the effects of externally applied force fields, one 
must not restrict the deterministic part of the Hamilto- 
nian to that of a free particle. For example, we will see 
that the presence of a constant magnetic field leads to a 
diffusive behavior in the sense that S| ~ t, a result with 
some relevance for magneto-transport theory. Secondly, 
the random part can be generalized to cover the case of a 
momentum-dependent (in other words, non-local) ran- 
dom potential, which is the continuous-system analog of 
off-diagonal disorder in lattice systems [|| . This is of in- 
terest, for example, to caricature the effective motion of 
a test particle due to inelastic scattering by the irregu- 
lar motion of other particles. It suggests itself also from 
the point of view of Hamiltonian mechanics. And lastly, 
we will show that the above-mentioned noise-induced re- 
sults are neither affected by quantum fluctuations nor do 
they depend on the initial state. 

In fact, it is the main purpose of the present Letter to 
demonstrate that there is a rather general class of Gaus- 
sian white-noise Hamiltonians for which one can obtain 
exact and explicit results on the averaged time evolution. 
Yet before we describe this class in detail, it seems ade- 
quate to comment on the representation we are going to 
use throughout. 

We consider a quantum-mechanical system which, for 
simplicity, has the Euclidean line R as its configuration 
space. The extension to the cZ-dimensional Euclidean 
space M. d is merely a matter of notation. Since the ran- 
dom part of the Hamiltonian will be allowed to depend 
on both position and momentum, it is convenient to char- 
acterize its properties in terms of those of an associated 
random function on classical phase space Ixi There- 
fore, and in order to treat the classical limit with low 
effort, it is only consequent to represent the quantum 
system entirely in phase space. The representation we 
choose is the one dating back to ideas of Weyl, Wigner, 
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and Moyal |Q, where a quantum operator / acting in 
the Hilbert space of square-integrable functions on M is 
represented uniquely by its symbol, that is, by the phase- 
space function 

f(p, q) := [ dr e^' h (q - r/2\ f \q + r/2) . (1) 



We recall that the symbol of the standardized commuta- 
tor (i/h)(fg — gf) of two operators / and g is the Moyal 
bracket of their symbols / and g 

[f,g](p,q) f(p,q) | sin |^ ( 9 p 9 q ~ d? 5 ^ j 9ip,q) ■ 

For non-polynomial / and g it is often advantageous to 
rewrite [/, g] as a Fourier-integral expression. Of course, 
in the classical limit, when Planck's constant ft tends to 
zero, the Moyal bracket reduces to the Poisson bracket. 
An observable, corresponding to a self-adjoint operator, 
has a real symbol a. A quantum state is represented by a 
Wigner density w which, by definition, is (27rft) _1 times 
the symbol of the corresponding "density matrix" , that 
is, of a positive unit-trace operator. The expectation 
value of the observable a in the state w is then given by 
the scalar product 



(w,a) 



dpdq w(p, q)a(p, q) 



We recall that \w(p,q)\ < (7rft) _1 and that (w,l) = 1. 
Moreover, (w,w) < (2ttH)~ 1 with equality if and only if 
w represents a pure quantum state. In the classical limit 
a quantum state w converges (weakly) to a probability 
density on phase space, that is, to a classical state. 

Now we are in a position to introduce a dynamics gov- 
erned by a Hamiltonian on phase space which we call 
white-noise Hamiltonian 



H{p, q) + N(p, q; t) 



(2) 



By definition, the Hamiltonian operator is obtained from 
(||) by inverting ([!]), that is, by Weyl ordering. For 
simplicity, the deterministic part H(j>, q) is supposed to 
be time-independent and at most quadratic in p and q. 
The random part N(p, q; t) is supposed to be a Gaussian 
white-noise field with mean zero and covariance 



N(p, q; t)N(p>, q>; t>) = C(p -p',q- q')S(t - t') . (3) 

Here the overbar denotes averaging with respect to the 
probability distribution of N and homogeneity is as- 
sumed just for brevity. By its probabilistic origin, the 
covariance function C may be any even phase-space func- 
tion with a non-negative (symplectic) Fourier transform 



C{x,k) 



dpdq 



For later purpose we will assume that C is sufficient- 
ly smooth, equivalently, that the probability density 
C(x, fc)/C(0, 0) has moments of sufficiently high order. 

In the Schrodinger picture the time evolution wq ^ w t 
of a given initial state wq is determined by the stochastic 
quantum Liouville equation || associated with 



d t Wt = [wt,H} + [w t ,N(t)} 



(4) 



Here the bracket [wt , H] is in fact a Poisson bracket due 
to the quadratic nature of H. In order to derive an equa- 
tion of motion for the averaged state wi(p, q) ■= wt(p, q) 
from (Q), we follow essentially the earlier treatments in 
HH and perform a functional integration by parts with 
respect to the Gaussian average M. In doing so, we 
think of the Dirac delta function in (|3|) as being approxi- 
mated by a sequence of smooth covariance functions with 
correlation time tending to zero, which amounts to the 
Stratonovich interpretation of (||). The final result 
can be cast into the form of a linear integro-diffcrcntial 
equation 

d t uh(p, q) = [wi, H] (p, q) 

+ -rj / dxdk C(x, k) {wl[p + Hk,q + Hx) - wi(p, q)} 

(5) 

which is valid for t > and has to be supplemented by 
the initial condition Wq = wq. Several remarks are in 
order: 

• Eq. (|^) is a substantial generalization of the main re- 
sult of H|. In the special case of a free deterministic part, 
H = p 2 /2m, and a momentum-independent random po- 
tential, C(x, k) — S(x)j(k), it reduces to an equation 
which is equivalent to Eq. (8) in ||. Furthermore, the 
subsequent treatment of Eq. (8) in || is restricted to a 
pure initial state represented by a joint Gaussian wq with 

(Tq ^ 0, and 



0, (w ,q 2 ) 



C(p, q) cos(xp — kq) >0 



(w ,p) = {w 0l q} = {w ,pq) 
(w ,p 2 ) = (ft/2a ) 2 . _ 

• The averaged time evolution T t : wq >— > w t given by 
(||) provides an example of a quantum-dynamical semi- 
group 1 1 1 1^] which is monotone mixing-increasing. By 
this we mean that T t maps Wigner densities linearly 
to Wigner densities, T t o %i = T t +t' for all t, t' > 0, 
Tq = identity, and d t (7~ t (wo) ,Tt(wo)) < (with equal- 
ity only in the uninteresting case where the covariance 
function C is a constant). The inequality follows from 
scalar multiplication of (|5|) by Wt = T t (wo), by observ- 
ing (Wt, \wt, H]) = 0, the Cauchy-Schwarz inequality and 
C > 0. To summarize, the average over randomness has 
turned the fully reversible quantum Liouville equation 
(^) into the equation (^J) with coherence-destructing ir- 
reversible behavior. 

• Interestingly enough, modified quantum-dynamical 
equations similar to (^) are discussed in very different 
branches of physics. These include quantum theories not 
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only of certain disordered systems, but also of the mea- 
surement process , of Markovian transport , and 
of the evaporation of black holes H] . 

• As for the generator of the semigroup, one deduces 
from Eq. || that 

% = e -*C*+*0 ( 6 ) 
with the unperturbed Liouville operator 

C := (d p H)d q - (d q H)d p (7) 
and the noise-induced, irreversibility causing operator 
M := h- 2 {C{0, 0) - C(-ihd q , ihd p )} . 

• For the derivation of explicit results it is often useful 
to isolate the unperturbed time evolution in (|^) according 
to standard perturbation theory 



% = e~ tc exp \ - J ds e sC N e" s£ 



(8) 



We note that due to the quadratic nature of H the oper- 
ators e sC Af e~ sC commute at different times s and can 
be written more explicitly as 

e sc e -sc = h -2{ C (0, o) _ c(-ih)C s , ihX s )} . (9) 

Here we have introduced the time-dependent first-order 
differential operators IC S := (d q e~ sC q)d q + (d q e~~ sC p)d p 
and X s := (d p e~ sC p)d p + (d p e~ sC q)d q , whose coefficients 
are obtained from the phase-space trajectory of the un- 
perturbed problem as indicated, and do not depend on p 
and q. 

• By using (||) and (||) it is straightforward to de- 
rive a Fourier-integral expression for the integral kernel 
Tt(p, q\p' , q') of %, which is the solution of (0) with ini- 
tial condition Tq(p, q\p', q') = 5(p — p')5(q — q'). Since we 
will not need this expression below, we omit it. 

• The quantum-dynamical semigroup T t admits also 
a purely classical interpretation. This is because it is 
positivity preserving, which follows from (||) and the fact 
that C(—ihd q ,ihd p ) and e ±sC are positivity preserving. 
Therefore, % maps classical states to classical states. 
In other words, its integral kernel T t (p, q\p' , q') can be 
interpreted as the transition density of a stationary, in 
general non-continuous, Markov process in phase space 
and Eq. @ may be viewed as the associated classical ki- 
netic equation. In fact, up to the drift arising from H, 
it is a linear Boltzmann equation with a (homogeneous) 
stochastic kernel It is used, e.g., in quasi-classical 
theories of charge transport in semiconductors . 

• Assuming an ^-independent covariance function C, 
one has the expansion 



where the three constants -Do. 2, -Di.i and Z?2,o a s defined 
through D^ v := (-id p )^(id q ) u c(o, 0)/fx\ v\ reflect the 
curvature of C at the origin and obey the inequalities 
D ,2 > and 4D 2 ,oA),2 > £>f,i due to C > 0. As a 
consequence, in the classical limit Eq. (||) reduces to a 
Fokker-Planck-type equation in phase space with drift 
and diffusion as given by ^ and (|l~0|). 

Now we return to the problem posed in the beginning 
of this Letter, namely to evaluate the averaged expecta- 
tion value {wt, a) of a simple observable a at time t, given 
the initial state Wo- For this purpose it is useful to switch 
to the Heisenberg picture according to 



(w t ,a) = (w u a) = (%(w ),a) =: (w ,T^(a)) . (11) 

The thus defined adjoint semigroup T t * can be obtained 
from (||) or (||) by reversing the sign of C. 

To be more specific, we first choose the deterministic 
part of the Hamiltonian (g) to be that of a free particle, 
H = p 2 /2m. Taking the observables p, q, p 2 , pq, q 2 , and 
q as examples, one then finds explicitly 





= P, 


V(q) = 


q + tp/m, 


(12a) 


V(p 2 ) 


= P 2 4 


- 2tD , 2 , 




(12b) 


V( P q) 


= p(q 


+ tp/m) - 


\-Wi A +t 2 D a a/m , 


(12c) 


V(q 2 ) 


= (q^ 


tp/m) 2 - 


^2Q 3 (t) , 


(12d) 


V(q 4 ) 


= (q^ 


tp/m) 4 - 


^12{(q + tp/m) 2 Q 3 (t) 





+2h 2 Q 5 (t) + (Q 3 (t)) 2 }. (12e) 



,fl ;1 CD,,-^,,-! m 1 v /u is a polyno- 



N- 



-D Qt2 d 2 - D^xdpdq - D 2fl 8 2 + 0(h 2 d 4 ) , (10) 



Here Q M (t) := ^Z. 
mial of (maximum) degree \i in time. Quantities such 
as the spatial variance of the averaged state E| := 
(wq, T t *(q 2 )) — (wo,T t *(q)) 2 or the averaged mean-square 
displacement A 2 :— (wo,T t *(q 2 ) — 2q%*{q) + q 2 ) at time 
t, may now be obtained immediately. 

The exact results (12) illustrate important features 
valid for general quadratic H: Observables which are 
linear in p and q are not affected by the white noise N. 
Moreover, noise-induced terms in averaged expectation 
values ([ll]) of quadratic observables are independent of 
the initial state. Assuming an ^.-independent covariance 
function C, noise-induced effects affected by quantum 
fluctuations occur only for observables of at least fourth 
order in p and q. However, very special situations are 
needed for quantum effects to show up in the leading 
term for long times. For example, taking the observables 
p n or (j", 4 < n integer, one must require the phase- 
space trajectories of H to grow exponentially in time. 
This remark contradicts certain expressions in since 
the underlying "correlation functions" considered there 
do not have a positive Fourier transform, and are there- 
fore physically insignificant. 

White-noise Hamiltonians may reveal a diffusive be- 
havior in a weak sense, that is, T t *(p 2 ) ~ t and/or 
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T t *(q 2 ) ~ t for long times. The simplest case for 
weak diffusion to occur in both momentum and po- 
sition corresponds to H = 0, as follows from (12t) 
and (12d) in the limit m — > oo. It occurs also if 
the deterministic part describes a harmonic oscillator 
and, more strikingly, in the case of a particle with elec- 
tric charge — e moving in the Euclidean plane M 2 un- 
der the influence of a perpendicular constant magnetic 
field of strength \mui/e\. In the latter case we choose 
H = {2m)- 1 {(pi - mujq 2 /2) 2 + (p 2 + m< J jq l /2) 2 } and, 
for the sake of brevity, we assume the covariance func- 
tion C(p,q) to depend only on the absolute values |p| 
and |q|. By a simple extension of the presented meth- 
ods to higher dimensions, one then finds for the averaged 
squared position at time t 



T t *(q 2 ) = (e t£ q) 1 



sin cut 



ajt 



1 



sin cut 

LUt 



d 2 qi C(0,0) 



(13) 



Here e t£ q is nothing but a cyclotron orbit associated with 
H. Li ke in lattice models ||, but unlike the free-particle 
case (12d), the leading term in (|l3j) for large t ;» uj^ 1 



is influenced by the noise in both momentum and posi- 
tion. Note also that the long-time behavior of 7^*(q 2 ), 
and hence of E 2 and A 2 , changes suddenly from a cubic 
to a linear behavior when turning on a magnetic field. 
When the fluctuating environment is modeled by a ran- 
dom part with a non-zero correlation time, we expect the 
same trend. However, the growth should be slower and, 
as a new feature, it should depend on the initial state. To 
verify these conjectures, non-perturbative methods are 
needed for controlling the long-time behavior of quanti- 
ties as E 2 and A 2 in these non-Markovian models. 

Not unexpectedly, except for very particular deter- 
ministic parts such as H — (p 2 — Dq 2 )/2m with D > 
-Do, 2/ ' D2fl 1 white-noise perturbations lead to a linear in- 
crease in time of the averaged energy 



t;(h) = h 



t (D 0>2 d*H 



DxsdpdqH 



D 2 , d 2 H) 



To compensate this effect and, if possible, to allow for 
an eventual approach to a stationary (equilibrium) state, 
dissipation has to be incorporated, typically by coupling 
the white-noise system to a heat bath in the spirit of jl7) . 
In the context of noise and dissipation, Accardi's program 
"Quantum Stochastic Mechanics" jl8| is very promising 
for a deeper understanding of quantum time evolutions. 
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